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Abstract 

Using the point fusion procedure we obtain the new integrable systems from the Elhptic 

■ Schlesinger system (ESS). These new systems have the pole orders higher than one in the 
I matrix of the Lax operator. Quadratic Poisson algebras on the phase space of the new 

■ systems generalize the Sklyanin algebras and have the graduated structure. 
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1 Introduction 



The examples of the systems with Lax operators having the pole orders higher than one were 
considered in papers [T] and [2]. These systems were obtained by point fusion procedure called 
I in algebra Inonu-Wigner ([3]) contraction. The main idea of this method consists of finding 

^ ' such decomposition of the variables, which gives us the pole order in the matrix of Lax operator 

being higher than one at some marked point. It implies the existence of Hamiltonian and 
Casimir numbers in this new system being the same as in the initial one. The main object of 
this procedure in the present paper is the Lax operator of the Elliptic Schlesinger system (ESS) 
^ ■ which was considered in [4j. 

! Historically ESS appeared as an approach to the decision of Riman's problem about searching 

the differential equations with regular singularities and intended monodromy data: Mj, j = 
Tjlj- ■ l,...n — > ^Mj), where ^ is the solution of initial linear system. In [S] L. Schlesinger had 

. considered the first order system of differential equation on CP^ 



S: (ft+E?.,^)* = o 

1 d,„1 = 

^ , and found the preserving conditions on matrix with respect to the changing the point 

^ I positions Xj. He obtained the system of first order differential equations for n > 3 matrices 

(j = 1, . . . , n), depending on n points Xk € CF^ 



O^S^ = , (k + j) , dk = a,, , (1.1) 

5,S'= = -^i^. (1.2) 

This system is non- autonomous Hamiltonian system and has the Hamiltonian form with respect 
to the linear (Lie- Poisson) brackets on sl(A^, C). The Hamiltonian 

]^k 
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defines the evolution with respect to the time Xk- 

For two by two matrices and four marked points the Schlesinger system is equivalent to the 
Painleve VI equation ([6],[1]). In this case the position of three points can be fixed as (0, 1, oo) 
while 2:4 play the role of an independent variable. 

If we replace CP^ by an elliptic curve, we define a similar system (the elliptic Schlesinger 
system (ESS)). In this case, in addition to the coordinates of the marked points, a new inde- 
pendent variable appears inevitably. It is the modular parameter of the curve, and thereby we 
have an additional new Hamiltonian. This system was introduced originally by Takasaki [7]. 
His derivation is based on the quasi-classical limit of the quantum SU{N) version of the XYZ 
model. In [4j was obtained as symplectic quotient of the symplectic space of connections of 
principle bundles of degree one over the elliptic curves with n marked points. This approach 
was previously developed in [8]. 

Let us give a short description of ESS. 

Elliptic Schlesinger System. Let us consider elliptic curve S,- = C/(Z -|- rZ) with the 
modular parameter r, (S(t) > 0) and 

— (2^1 1 • • • ) Xfi) ) Xj ^ xj^ , X}^ G Yjf 

be the divisor of non-coincident points with the condition 

^ Xj € (Z rZ) . (1.3) 

Consider the space Vn,N of n copies of the Lie coalgebra g* sl{N, C)*, related to the points 
of the divisor. 

where T" is the basis element of GL{N,C) (see Appendix B). 
Introduce operators acting from Vn,N to the dual (B^^iQj 

Ifcj : flfe 9j , S'^ ^ {Ikj)ySi^ , {Ikj)y = (py{xj - Xk) , (1.5) 
Jji : Q*j 9j , Sl^ ^ J^Si^ , = £^2(7) , (1-6) 
Jfej : 0fc 9j , S':^ ^ {Jkj)'ySif , {Jkjh = fyixj - Xk) (1.7) 

where ip^{x), -£2(7)5 f-y{x) are defined in Appendix B. 

The positions of the marked points xj G Dn and the modular parameter r are local coordi- 
nates in an open cell in the moduli space Aii^n of elliptic curves with n marked points. They 
play the role of times. 

Definition 1.1 . The elliptic Schlesinger system (ESS) is the consistent dynamical system on 
'Pn^N with independent variables from Aii^n 

djS'' = [IkjiS^), S"] , (k + j) , dk = d,, , (1.8) 



afcS'= = -^[I,fc(S^),S'=], (1.9) 

= E j,,(s'^)] + ^[s^ j,,(s^)] , (1.10) 



where the commutators are understand as the coadjoint action of Qj on q'j . 
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In the basis (a G zj^^) dEB the ESS takes the form 

dkSi= C{j,a)S';Si_^ip^{x,-Xk), (k^j), (1.11) 

45^= C(7,a)5^5L^5;^'99„_^(xfc-x,), (1.12) 

^-*^' = i E C(a,7) I J^5^^5^^(xfc-x,) + S:^S^^i?2(7)) • (1-13) 



Remark 1.1 . In the rational limit C^rriT ^ oo) U.ll\) and U.l^) pass to the standard 
Schlesinger system lll.l]) . (see \A.l^) ). 

As in the rational case the ESS has some fundamental properties The space "P^ ^ is Poisson 
with respect to the linear Lie-Poisson brackets on g* 

{SlSlh=5^^C{a,l3)Si^^ (1.14) 

ESS is a non-autonomous Hamiltonian system on Vn,N 

dkS^ = {Hk, }i , dk = d^, , (1, . . . ,n) , (1.15) 

drS^ = {Ho,S^}i, (1.16) 



where 



Hk = - Yi^kj{S')Sn) = - E E S'.Si^ip.ix, - xu) , (1.17) 



(1.18) 



1 

'27n 



Y E si^s'_^f,{xk - X,) +J2 E ^^^-7^2(7) 



The brackets (|1.14p are degenerate. The symplectic leaves are n copies of coadjoint orbits Oj 
(j = 1, . . . ,n) of SL(A^, C). Assume that all orbits are generic, and let c^{j) be corresponding 
Casimir functions of order /i (/i = 2, . . . , N). The phase space of ESS is 

■Rn,N ~ Vn,N/{c^{3) = c'^(i)o} -\{0, , (1.19) 

dim7e„,Ar = nA^(iV-l). (1.20) 

The ESS can be considered as a system of interacting non-autonomous SL(A^, C) Euler- Arnold 
tops, where operators (jl.Sp . (II. 6|) . (II. 7p play the role of the inverse inertia tensors. 
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2 Classical integrable systems obtained from ESS 

ESS in the case of two marked point and N = 2. The Lax operator of ESS in the case 
of n marked points has the following form: 

1 " 

Liz) = --E,iz)To + ^ Yl - • (2-21) 

where Tq and are basis elements of GL{N, C). 
Let us consider the modified Lax operator ([4j) 



n 



L3™-P(z) = 5oTo + S^Eiiz - x,)To + ^ 5>„(z - xj) , (2.22) 

j \ a / 

where we attribute to the marked points of the divisor Dn n copies of the GL(A^, C)-valued 
elements 



adding to this set a variable 5*0 € C. So, it defines 



P+^ = {5o, (5^,S^J = l,...,n)| ^5^ = 0}. 

i=i 

It is possible to obtain the equation of motion for ESS from the quadratic brackets on the 
space ^, extracting them from the classical exchange algebra 

{L3™«P(z), L9™«P(ii;)} = [r{z - u;), L9™"P(z) (g) L^^^^^^Cu;)] . (2.23) 

where r is the classical Belavin-Drinfeld r-matrix. 

Let us consider the first occurrence: N = 2 and n = 2 is a number of the marked points. 
Then the Lax operator takes the following form: 



L™(z) = [So + S^^Eiiz - xa) + S^Eiiz - Xk)) ao+ 

+ (^>"(^ - + - "^b)) c^a , (2.24) 

a 

+ S', = 0. 

Note, that for N = 2 the basis is proportional to the basis of the Pauli matrices. The 
dimension of the phase space i?2,2 is four. The Hamiltonians of this system are Sq and Sq. The 
quadratic brackets: 

drS-^ = ^{S-^,So}= (2.25) 



(^2(7) - E2{(3)) - ieaf3yy^^{Xab)SpS^ + ieaf3^^f3{Xab)SpS: 
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= = ^ {5^, S^} = (2.26) 
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System obtained via two point fusion. Let us fulfil the following coordinate transformation 
and decomposition of the variables S ([Ij): 

Xb = Xa + e, Sq = c^ qSq + cI^ iSqE = c^^^Sq + cI^^^SqE , (2.27) 

where c are some coefficients. Taking the limit e — > and putting some additional conditions 
on coefficients (the requirements of the absence of singularities) we get the new Lax operator 

-^fusion i^) 

Lfusioniz) = {So + S^oE2{z)) ^tq + {s^aVa{z) + Sl^'^{z)^ , (2.28) 
where we put the position Xa equal to 0. The dimension of the phase space 

'Pf,2 = {Sq,sI,s^,sD 

is 8 and it is equal to the dimension of ^ before point fusion. 

To find the exchange relations between the new variables S we consider the equation (12.230 
and the new Lax operator (j2.28p . These exchange relations are the coefficients at the function 
products 1, E2{z), (pa{z), fa{z) and 1, E2{w), (pa{w), ifaiw)- So we get 

Proposition 2.1 . The space 7^/2,2 "is Poisson with respect to the quadratic brackets 

{5o,5o'}=0, (2.29) 



^ai f ~ '^i^ap-yS^S^ — 2zeQ/3^£'2(7)S'oS'° 



(2.30) 



1 cl 

7 ' 



(2.31) 



[Si, So] = 2ie^p^ (^2(7) - E2{P)) SoS^ + 2ie„^^^2(a) (^2(7) - ^2(/3)) S^S 
{Si So} = 2ie^^^E2{P)S''pS] - 2ie^^^E2{j)S}S^ , 

{S^ Si} = 2ie^^, (i?2(7) - E2m S}S}^ , 
{S^i So} = ~2i£af3^S^S}^ + 2ieaf)^S^S^ . 

Proof. 

Let us consider the classical Poisson brackets 

{L{z),L{w)} = [r{z - w),L{z) O L{w)] , (2.32) 
where the Lax operator has the following form: 

L{z) = Lfusioniz) = {So + SlE2{z)) ao + Y, {sl^c.{z) + S^(^l(z)) , (2.33) 



In the l.h.s. for the matrix element (Tq ® CTjj in (IA.24P we have the sum consisted of the 
following terms: 

{'S'a; Sl}Lpa{z)ipp{w), {5"°, Sl]Lpa{z)Lp p{w) , 
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and for the matrix element (Tq /: 

{SlSo}^a{z), {SlS'o}Va{z)E2{w), {5^, 5o}^l(z), {Sl S^,}cp'^{z)E2{w). 
In the r.h.s. we have for the same matrix elements the sums consisted of the following terms: 





2iea3jS^So [(pa[Z - W)ip^{w) 


- "Pisiz - 


■w)(p^{z)) , 






2l£ap^b'^bQ [(fia[Z - W)(fij{w)E2{z) 


f 

- ^p{z - 


w)(pj{z)E2{w)) , 






ZlEa/S-yJ^OO [^a{Z - W)ip^[w) 


- ipp{z - 


■ w)<p^[z) j , 






2i£aPjS}^SQ (^(pa{z - w)(p'^{w)E2{z) 


- ^niz - 


- w)(p'^{z)E2{wfj 




and 












2i£a^^S^S'^ {ip^{z - w)ip^{z)ip^{w) ■ 


- ff^{z - 


w)(p^{z)(P(^{w)) , 






2ieap^S}^Sp (^(p^{z - w)(pi3{z)(p^{w) ■ 


-(Pisiz- 


w)(p'^iz)ipi3{wfj , 






218^(3-^3^81 (^(p.y{z - w)ip'p{z)ip^{w) ■ 


- 'Pf^iz- 


w)(p^{z)ip'p{w)^ , 






2ieapjS}^S'^ (^ip^{z - w)(pp{z)(p^(w) 


-(p/3iz- 


-w)^p^{z)^pp{w)) 




As a result of the expansion in function product for the matrix element ctq (8) 


we get: 



2i£a(j^S^So {ipa{z - w)(p^{w) - ipp{z - w)ipj{z)) = 2i£af3^(pa{z)ipp{w)S^So, 
2i£al3-iS^Sl {(Pa{z - w)(p^{w)E2{z) - ip/3{z - w)ipj{z) E2{w)) = 

= 2i£ai3^ (^(pa{z)(pp{w)E2i'y) - (p'^{z)(p'i3{wfj S^Sl, 

2i£al3^S]SQ {^^a{z - W)ip^{w) - ip/3{z - w)ip'^{zj^ = 
= 2i£apj (^^p'^{z)(pp{w) + ipaiz)ip'p{wj^ S^Sq, 

2ieaP'yS^SQ (^ipaiz - w)(p^{w)E2{z) - ipp{z - w)(p^{z)E2{wj^ = 

= 2iea/37 [v'aiz)^(3{w)E2{P) + ipa{z)ip'i3{w)E2{aj^ S}^Sq. 
and for the matrix element Ga® I'- 

2ieaf3'yS^S'^ {(pj{z - w)ipp{z)ipj{w) - ipp{z - w)ipj{z)ipf3{w)) = 
= 2ie^p^^a{z) (i?2(7) - E2m S^^Sl 

2ieai3yS}^S'^p (ipy{z - w)ipfs{z)ip'^{w) - ipp{z - w)ip^{z)Lpp{w) \ = 
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liEafj^S^^Sl (^ip^{z - w)Lpp{z)ip^{w) - iffsiz - w)cpy{z)ip'p{w)^ = 

2iea/3^S}^Sl (^if^iz - w)ip'p{z)ip'^{w) - Lpp{z - w)ip^{z)Lpp{w)^ = 
= 2ie^p^ {-^a{z)E2{a) (^2(/3) - ^2(7)) + Mz)E2{w) (^2(/3) - £^2(7))) S},S}. 

After regrouping the terms we get (j2.29p , ()2.30p and (j2.3ip as the coefficients at the function 
products 1, E2iz), ifaiz), ^a{z) and 1, E2{w), ifaiw), ip'a{w). □ 

Let us write these exchange relations in the form of the tables (1) and (2): 



Table 1. Exchange relations Sa, Sp 



2i£f^f^^- 






ci cO 


cl cl 




+1 















+1 












+1 















+ 1 






and 

Table 2. Exchange relations Sa, Sq 



2i£Q,jjj- 




cO cl 


cl cO 


cl cl 


{Sa, So} 












{5*0, 5*0} 























+J7/3 







-1 


+1 
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where J a = -£'2(0), J^p = £"2(7) — -E'2(/3)- The Jacobi identity for Vf^2 follows from the 
classical Yang-Baxter equation for r matrix. Note that the Poisson algebra Vf^2 come to the 
Sklyanin algebra ([9j), if we put Sq =0, 5^ = 0. 

The quadratic brackets are not degenerate on the orbits. To describe the system we define 
the Casimir functions and Hamiltonians. The equation for the spectral curve has the following 
form: 

det{L{z) - XI) = 0, 

- Tr{L{z))\ + detL{z) = 0, 

where det is determinant. The coefficients TrL{z) and detL{z) of this equation define the 
Casimir functions and Hamiltonians of the system. TrL{z) and detL{z) are doubly periodic 
functions and they can be decomposed into the basis of Eisenstien functions. 

]^TrL{z) = So + E2{z)Sl, 

detL{z) = Co + CiEi{z) + C2^2(^) + C?,E'^{z) + C^E^iz). 
So, we have two Hamiltonians S'o, 5q and four the Casimir functions (Ci = 0): 

Co = Sl-Ar,l{Slf + E2{a){Slf+ 
+ [{E2{a)f - {E,{a)fE2{a)+E,{a)E'^{a) + ^4' (a)) (5^)^, 

C2 = S,Sl + 4r?i(5i)2 - {Sy + {{Er{a)f - E2{a)) {Sl,f, 

C3 = —SaS^, 

c. = \ {{slf - {slf) . 

To calculate the Casimir functions we use (IA.7P - (IA.10|1 of Appendix A. The dimension of 
the phase space TZf^2 (symplectic leaves) is equal to 4. 

nf^2-Vf,2/{Ci,i = T;4} 

In terms of the quadratic brackets the equations of motion have the following form: 

dto^a = 2^311,30} = 
= isaf^^ {E2{l) - E2m + ie^f,^E2{a) (i?2(7) - i?2(/?)) S^S}^ , 

(2.34) 

dt.S^a = US'a, s'o} = iSo^M (^2(7) - E2m Sp}^ , 
dt^S^ = \{S\, Sq} = —iEap-yS^S}^ + i£al3jS^S!^ . 
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Graduation and the systems via three point fusion. Let us do the following observation. 
One can see from the tables (3) and (4) that coefficients at SS products have the certain 

dependence of the same products. Consider the bracket j'S'lJj'S'J^I, we note that SqS^ differs 

from SqS^ by one of the multipliers. The coefficients take the values +1 and — J-y correspondingly. 

For the bracket |s'J^,5^| we have the same situation: the product SqS}^ differs from SqS}^ one 

by the same multipliers Sq and Sq. The coefficients take the values +1 and —J/3. Now the 
coefficient at 5q5^ depend on /3. Taking into account Ju = E2{u) = —-^Ei{u), we see that 
each function J has the degree of the operator ^ proportional to 1, but it takes the different 
value. It is possible to say that the function +1 has the degree of the operator ^ equal to 0. 
It is possible to see the regularization in the coefficient positions for the other brackets. So, we 
can define the notion of the graduation in the following manner. We consider the Lax operator 
Lfusion{z) (I2.28|) and assign that the graduation of variables S - grad{S) = [S] is proportional 
to the modulus of the pole order of the function at S in (j2.28p . It is possible to assume that 
in the expression of the each Poisson bracket all the terms have the equal graduation. Here we 
take into account that the graduation of the coefficients (which is proportional to the degree of 
the operator) is opposite in sign to the SS products. The graduation of each term is the sum 
of the graduations: 

[7^515!] = [J^] + [5i] + [5i]. 

From this approach the interesting fact is that the operation of taking Poisson bracket acquire 
the graduation too. 

[{5°,5°}] = [{,}] + [5°] + [5°]. 

Now let us compose the linear equations. In a simplest case of the Sclyanin algebra ([9]), if we 
put 5^ and Sq equal to 0, we get the following equations: 

2b+[{,}]=a + b, 
a + b+[{,}]=d + 2b, ^^-""^^ 

where we put [Sq] = a, [S^] = b, [J^] = d. We get the following important relation 2[{, }] = d. 
Note that we get it without using the value of graduations a and b. 

Let us consider now the three point fusion. In this case the Lax operator has the following 
form: 

= (So + S'qE2{z) + SiE',{z)) ao + Y, {s>a{z) + + Sl^l{z)) . (2.36) 

a 

Analyzing tables (3) and (4) it is possible to conclude that the "structural blocks" for the 
coefficients at SS are the product of the function J. Using the graduation we can write the 
presumable form of the exchange relations. Let us consider for example {S^, Sq}: 

{'S'a; Sq] = 2i£af3-y [JjjjSoS^ + JaJjjjSjjS}^ + C2fiS'^S^ + Co^2S^S'^ + C2,2S'fjS'^) . (2.37) 

Put [So] = 0, [5°] = 1, = 2, [Si] = 2, [S^] = 3, [S^] = 3 in accordance with notation above, 
we get [Ja] = —2 and [{S^, Sq]] = 0. So there are all allowed terms in (12.370 . For example it is 
not possible to write the term with S^S'^. Its graduation is equal to 5, and we cannot pick out 
the coefficient consisted of J-function product because the graduation of the whole term is not 
equal to for each case. Compared the results of the explicit calculation ((Table 3), (Table 4)), 
we get the coefficients. The graduation shows the possible positions of the coefficients. Note 
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that the coefficients at SS are the invariants of the transformations which do not change the 
Poisson brackets. 

Table 3. Exchange relations Sa, 





5*0 5*^ 




SqS^ 






q2 cO 


SqS^ 




C2 c2 


/ cO cO \ 


+1 





-J^ 


















fcl c-ol 





+1 

























+1 





-Ji 





-./- ( 


















+1 





+1 





-~{.Ja + Jl3) 








/ CO 02 \ 














+1 












f C2 qO \ 














+1 












fcl c2\ 

















+1 








— Jal3 


( c2 cl\ 

















+1 








Jal3 


/ C2 02 \ 




















1 

2 


+h 






Table 4. Exchange relations Sq,, 





CO CO 


SffS}^ 


SpS^ 


cl cl 


cO c2 


C2 cO 


cl C2 


C2 cl 


q2 c2 
















JaJ'y/3 








(J/3 J7 - Ja)Jj0 


{'S'i, So} 




















J 8 Jya 





{S'°, /Sq} 






























-1 


+1 


















{s*^, 5*0} 














































1 T 

2 "^7/3 


1 T 

2^7/3 



















-1 


+1 
























-1 


+1 










{5^, S'o} 




















1 

2 








As the result we get the following 
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Proposition 2.2 . The space Vj^^ is Poisson with respect to the corresponding quadratic brack- 
ets Table. 3 and Table. 4. 



3 Conclusion. 

Although the graduation clarifies the algebraic structure of the phase space by indicating the 
possible positions of the coefficients at SS product but it is not possible to write the explicit 
forms of these coefficients. The other question is about the representation of the coefficients by 
J-function product. The goal of the further researching can be simulate the structure of the 
phase space without using the calculation and abstracting from the elliptic dependence. 

Acknowledgments. Author would like to thank A. Levin for fruitful discussions and M. 
Olshanctsky for the suggested theme, important discussion and attention during the writing 
of article. The work was partly supported by grants RFBR-06-02- 17381, NSch-8065-2006.2, 
RFBR-06-01-92054-KE and by Federal Nuclear Energy Agency. 

4 Appendix. 

4.1 Appendix A. Elliptic functions 

We assume that q = exp27riT, where r is the modular parameter of the elliptic curve Et-. 
The basic element is the theta function: 

■d{z\T) = q^ ^{-l)'^e{^n{n + 1)t + nz) = (e = exp27rz) (A.l) 
nez 

The Eisenstein functions 

Ei{z\t) = d^\og^{z\T), Ei{z\T)r^^-2riiz, (A.2) 

where 

''^^"^ = ^^' ^(-) = 11(1 - '?'^) • (A.3) 

' n>0 

is the Dedekind function. 

E2{z\t) = -d,Ei{z\T) = d'^^log^{z\T), E2{z\t) ~ ^ + 2r?i . (A.4) 

Relation to the Weierstrass functions 

Ciz,T) = Eiiz,T) + 2r]i{T)z, piz,T)=E2{z,T)-2rjiiT). (A.5) 
The highest Eisenstein functions 

(P{u, z)(P{-u, z) = E2{z) - E2{u) , (A.7) 
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(j)iu,z)(l)'{-u,z) = + (^Ei{u)E2{u) + ^E2{u)'^ - Ei{u)E2{z) + ^^^(z) 

(t^{u,z)(t){-u,z) = - {ei{u)E2{u) + ^^2(n)') + Ei{u)e'^{z) + 

4> {u,z)(t) {-u,z) = 
-{E2{u)f + {Ei{u)fE2{u) + Ei{u)E2{u) + ^E2{u) 

+ {-{Ei{u)f + E2{u)) E2{z) + ^-E2{z) . 
The next important function is 

^, , i9(u + z)t?'(0) 

(Piu, Z) = ^ s • 

1J[U)V[Z) 

4){u, z) = 4>{z, u) , 4>{—u, —z) = —(j){u, z) . 

It has a pole at z = and 

ct>{u, z) = ^ + Ei{u) + ^{El{u) - p{u)) + ... . 

du4>{u,z) = 4){u,z){Ei{u + z) - Ei{u)) . 
dz(l)iu,z) = (j){u,z){Ei{u + z) - Ei{z)) . 



Heat equation 
Quasi-periodicity 



\im \n du4>{u, z) = —E2{u). 



dr(t>{u,w) - ■^dudw(f>{u,w) = 0. 



^{z + 1) = -^{z) , i9{z + r) = -q-2e-^^''^z) , 
Ei{z + 1) = Ei{z) , Ei{z + r) = Ei{z) - 2m , 
E2{z + l) = E2{z), E2{z + t) = E2{z) , 
(j){u, z + l) = 4>{u, z) , z + t) = e~^™0(n, z) . 

dufpiu, z + l) = ducpiu, z) , du(t){u, z + t) = e"^™ z) - 2TTut){u, z) . 
The Fay three-section formula: 

(p{ui,Zi)(p{u2, Z2) - 4){UI + U2, Zi)(j){u2,Z2 - Zi) - (j){ui + U2, Z2)(piui, Zi - Z2) = . 

From IIA.I^} and \A.23\} we have: 
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(I){ui,z)(p{u2,z) = (l){ui + U2, z){Ei{ui) + Ei{u2) - Ei{ui + U2 + z) + Ei{z)) . (A.24) 

Particular cases of this formula are the functional equations 

4>{u, z)dy(f){v, z) - (t>{v, z)du(t>{u, z) = {E2{v) - E2{u))(t){u + v,z) , (A.25) 

0(n,zi)(/.(-n,Z2) = 0(n,zi - Z2){-Ei{zi) + Ei{z2) - Ei{u) + Ei{u + zi - Z2)) = (A.26) 
= 0(n, zi - Z2){-Ei{zi) + Ei{z2) + du(t){u,Z2 - zi)) , 

0(n, z)(t){-u, z) = E2{z) - E2{u) . (A.27) 



z — w)4>{ui — V, z)<j){u2 + v,w) — (j){ui — U2 — V, z — w)4>{u2 + V, z)(p{ui — V , w) = (A. 28) 

(j){ui,z)(j){u2,w)f{ui,U2,v) , 

where 

f(ni,U2,u) = Ei{v) - Ei{ui - U2-v)+Ei{ui - v) - Ei{u2 + v) . (A. 29) 

One can rewrite the last function as 

^'{0MuiMu2Mu2-ui + 2v) 
t(ui,U2,v)- ^(^^ _ + - ni + t;)i?(z;) • ^ ' 

Using ()A.2p . ()A.4p . (|A.13p one can derive from ()A.28p some important particular cases. One 
of them corresponding to v = ui (or v = —U2), is the Fay identity (lA.23p . Another particular 
case comes from ui = (or U2 = u): 

(j){v, z — 'w)(j){—v, z)(f){u + v^w) — (j){—u — v,z — w)(i){u + V, z)(j){—v, w) = (A. 31) 

(l){ui,z){E2{u + v) - E2(v)) . 
If n2 —V then (IA.28h in the first non-trivial order take the form for ui = a, U2 = (3 

4>{-l3, z - w)Ei{w)4>{a + /?, z) - 4>{a, z - w)Ei{z)(j){a + /3,w) = (A.32) 

(/.(a, z)<P{/3, w){Ei{a) + Ei{/3) -Ei{a + (3)). 

4.2 Appendix B. Lie algebra sl(A^, C), Group GL{J\f,C) and elliptic functions 

Introduce the notation 

27ri 

eN{z) = exp(— z) 

and two matrices 

Q = diag(eAr(l), . . . , eN{m), . . . , 1) (B.l) 

A = , {j = 1,...,N , modN). (B.2) 

Let 

Zj^^ = (Z/iVZ e Z/NZ) , zj^^) = zj^^ \ (0, 0) (B.3) 

be the two-dimensional lattice of order A^^ and A^^ — 1 correspondingly. The matrices Q^^A"^, 

(2) ~ (2) 

a = {ai,a2) G Z]y generate a basis in the group GL(A^, C), while Q'^^A"'^, a = (01,02) € Z]y 
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generate a basis in the Lie algebra sl(A^, C). More exactly, we introduce the following basis in 

(2) 

GL(A^, C). Consider the projective representation of Z)^ in GL(A^, C) 

a^T, = -^e^(^)Q-iA«^ , (B.4) 
Ztti 2 

TaTb = ^eN{-^^)Ta+b, (0x6 = 0162-0261). (B.5) 

^^^^ ^ non-trivial two-cocycle in Z2Ar). The matrices Tq,, a € 

generate a basis in sl(A^, C). It follows from (|B.5p that 

[T^,Tp] = C{a,P)To,+p, (B.6) 

where C(a,/5) = ■^sin-^(a x /3) are the structure constants of sl(A^, C). 
For N = 2 the basis Tq, is proportional to the basis of the Pauli matrices: 

^(1,0) = — ^^3 , 7^(0 1) = — CTl , T(i 1) = — CT2 . 

The Lie coalgebra g* = sl(A^, C) has the dual basis 

fl* = {S = ^^ = ^^-7> {T^tP) = 5-J. (B.7) 

It follows from ()B.6p that g* is a Poisson space with the linear brackets 

{So.,Sp} = C{a,P)S^+p. (B.8) 
The coadjoint action in these basises takes the form 

ad^J^ = C(a,/3)r+^. (B.9) 
Let 7 = '^^"ly^"^ - Then introduce the following constants on Z^^^: 

^(7) = ^(^^^^) , E,{^) = E,C-^±^) , E,{^) = E,C-^±^) , (B.IO) 

</.^(z) = </)(7,z), (B.ll) 
if^{z) = eN{^2z)(t)^{z) , (B.12) 

"PiA^) = ^N{l2z)(p{ri + ^^^'^ '^) • 
They have the following quasi-periodicities 

(f^{z + 1) = eN{-^2)f-i{z) , (f^{z + t) = eN{--ii)ip^{z) , (B.14) 

+ 1) = eAr(72)(/7^,^(z) , (/j^,^ (z + r) = Gat (-71 - r/) (/j^,^ (z) , (B.15) 
The important formulas with i^aiz) 

= -^^{z)^p{z) = ip^iz){Ei{-f + z)- Ei{z) - i?i(7)) , (B.16) 
ip''{z) = ip^{z){E[{a) + E[m - 2ip^{z)E[{z) . (B.17) 
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